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Abstract 

We analyze, by means of Mclnikov method, the possibihty of modifying the threshold of homoclinic chaos 
in general 1-dimensional problems, by introducing small periodic resonant modulations. We indicate in 
particular a prescription in order to increase the threshold (i.e. to prevent chaos), and consider then its 
application to the bistable Dufhng-Holmes potential. All results are confirmed both by numerical and by 
analog simulations, showing that small modulations can in fact sensibly influence the onset of chaos. 
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The problem of "controlling chaos" has received great attention in recent years [1-4]. Even if, for 
many authors, this term means in general the stabilization of unstable orbits, we are concerned here more 
specifically with the modification of the threshold for the onset of chaos. One of the simplest and most 
interesting methods used in order to modify (possibly to increase) the threshold of chaos which appears in the 
presence of homoclinic orbits, is the introduction of a periodic modulation of the parameters describing the 
unperturbed potential. This possibility has been analyzed for both the Duffing-Holmes and the Josephson- 
junction potential [2,3]. 

In the present note, we want to propose some generalizations of this idea: we will consider first a 
generic 1-dimensional problem by means of Melnikov theory [5-6], and obtain a simple criterion about the 
"correct" choice of this modulation; next, we will apply this indication to the bistable Dufiing-Holmes 
potential, introducing two independent modulations: this will allow us to confirm both the general results 
and the possibility of modifying the onset of chaos. All the theoretical discussion is well supported both by 
numerical and by analog simulations. 

Let us start by considering the general case of a 1-dimensional " equation of motion" for the real variable 
X = x{t): 



where for e = we have the standard problem of a periodically forced and linearly damped motion, whereas 
the additional term with e ^ takes into account the presence of general modulating terms. We assume 
as usual / = —dV/dx, where the unperturbed potential V = V{x) has a maximum point a,t x = xq and a 
homoclinic orbit, that we indicate by 



doubly asymptotic iov t ±oo to xq. In order to obtain a theoretical estimation of the effect of the last 
term of Eq. (1) to the threshold of chaos, let us write down the Melnikov function M{to) for the problem 
(1): taking into account that q{t) = q{—t), it is easily seen that M{to) acquires the form 



X = f{x) — 6x + 'ycoscot + eg{x) cos{flt + 9) 



(1) 



q = q{t) 



(2) 




(3) 



/' + 00 



esini 



.(l^to + ^) / q{t) g{fl{t)) sixiO^t dt = 6 Js -\- 'yJ-ySintuto -\- eJcSin{Qto -\- 0) 



J — oo 



2 



At the resonant case, i.e. when uj = Ct, this can be written 



-M{to) = SJs + jK sm{(jto + a) (4) 



with 



K 



1/2 



The condition for the appearance of chaos according to Melnikov criterion, i.e. that M{to) has simple zeros, 
becomes now (it is not restrictive to assume 7 > 0, whereas ^ > for physical reasons, and clearly Js > 0): 

> SJs (5) 

With fixed damping 5, we can then use the modulation term in order to modify the threshold of chaos: here, 
this amounts to modify the range of the forcing amplitudes 7 which do not produce chaotic responses. Now, 
according to Eqs. (4-5), the maximum increasing of this threshold is obtained by choosing 6 = and the 
sign of e according to the following prescription: 

e > (resp. < 0) if J-yJ^ < (resp. > 0) (6) 

(or, which is the same, e > in any case, ^ = if J-y Je < and 9 = it if J-yJ^ > 0). 

With this choice for e, and observing that the amplitude |e| of the modulation is usually very small (here 
we only require |e| < 7|J-y/Je|), then the above condition (5) becomes 

J, 



SJs ^ 



(7) 



or also 



where 



is the "Melnikov ratio" in the case e = of no modulation. 

Then Melnikov theory predicts that the presence of modulation produces, if the phase of modulation is 
chosen according to the above rule, an increasing (proportional to the quantity \Je/J^\) of the threshold of 
chaos. 
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Before further discussing this result, let us remark incidentally that in the case of non-resonant modula- 
tions, u) ^Q, Eq.(3) suggests that one may expect (at most after some time delay of the order ~ l/(w — il)) 
a "in-phasc" contribution of both forcing and modulating terms ~ 7I J-^j + |eJe|. Then, the above arguments 
indicate that a lowering of the threshold, favouring the onset of chaos, is to be expected in this case. A 
careful analysis of nonresonant modulations in the case of DufRng- Holmes potential can be found in ref.[2]. 

We want now to provide a precise test of the above results by checking them in the case of Duffing-Holmes 
potential. We find it convenient to introduce a small generalization of the above discussion, by choosing 
the perturbation g{x) in (1) in the form of two independent terms modulating both the linear and the cubic 
components of the force: this requires the presence of two modulation parameters €1,63, as made in [4], 
where a DufBng-Holmes - type potential was obtained by means of a magnetoelastic equipment. Precisely, 
we consider the equation 

X = —Ax^ (1 + £3 cos fist) + Bx{l + ei cos flit) — 5x + ^ cos u>t (9) 

where A,B,'y,6 > 0, and the signs of 61,63 are for the moment undetermined. As well known, the unper- 
turbed DufBng-Holmes potential possesses two homoclinic orbits, given by 

gW=±(^— j sechi\/B (10) 

where the signs -|- and — denote respectively the orbit surrounding the minimum of the potential at x^~^^ = 
-V^BjA and at x^~^ = —^BjA. All integrals appearing in Melnikov function can be evaluated exactly; 

in particular, the integrals ji^"* and ji^ take the same value when evaluated respectively along the orbit 
q^'^^it) and along q^~^{i), the same is true for and ji^'': we have 



= -2—— / sinh t/B sech^ tVB sin flit dt = cosech^i < (11) 

A A 2VB 

J^(±) = r sinh t\/B sech^ t\/B sinilst dt =^{fll+ 4B) cosech^!^ > (11') 

A J_oo QA ly B 



(notice that the results in Eq. (9) of ref.[2] and Eq. (4) of ref. [4] are uncorrect). One has also J^^^ = ^' 
whereas J^^-* = — J^"-*: 



Thus, Melnikov condition for the appearance of chaos, in the resonant case w = Qi = O3, can be finally 
written in the form 

J>i?(±)(ei,e3)=i?oTeiri±e3r3 (13) 
where again ± distinguish between the two homoclinic orbits q^'^\t), and 

Ro = ^ cosh — ■= (14) 

is the well known ratio for the unperturbed Duflang-Holmes potential [6], and 

n = — ^cotgh — -j= , rz = 1= — cotgh — -= (14') 

Therefore, we can conclude: If the motion occurs near the homoclinic orbit q'^^\t), the best choice in 
order to prevent chaos (by increasing the threshold) is 

ei < and €3 > (15) 

This choice however favours the onset of chaos when the motion is in the potential well around = 



— ^JB/A. The opposite choice ei > 0, es < would produce exactly opposite effects. All these results agree 
with our above discussion (cf. the signs of J^^' and Jef\jet^)- 
Numerically, with e.g. A = B = w = 1, 6 = 0.25, we obtain 

n = 3.08 , rs = 2.57 (16) 

The relatively large numerical values of these coefficients (compared with Rq = 0.753) show that the role of 
modulation in the Duffing-Holmes potential is really important: according to Eqs.(13,16), one may in fact 
expect that very small £1,63 may considerably influence the onset of chaos. Another interesting remark is 
that the introduction of the modulation (clearly the effect is present also choosing one of the two ei, £3 equal 
to zero) produces a sort of " dynamical asymmetry" between the two potential wells. Let us emphasize that 
it had been already remarked [7] that a small "geometrical" asymmetry in the double- well potential sensibly 
modifies the thresholds of chaos in the two wells. Precisely, considering an asymmetric potential 

V{x) = ^x^ - |a;2 + px (17) 



it can be shown that the Melnikov ratios i?^^^(/3), at the first order in the asymmetry parameter /J, are given 



by 



R^^\/3) = Ro ± /3p 



(18) 



where 




(19) 



With A = B= uj = las before, one has p = 1.35, showing that the two effects are in fact comparable. 

We have tested the above theoretical results for the Duffing- Holmes potential by means of both numerical 
and analog simulations. The agreement is globally good enough. For instance, choosing A = B = ui = 1, 
and S = 0.2, 7 = 0.2 (i.e. largely within the chaotic region if no modulation would be present, cf. [6]), it is 
sufficient to insert a modulation with ei = —0.05, £3 = +0.04 in order to obtain a periodic response. Fig. 
1 shows some of these periodic motions, oscillating around a;^+^ = 1, obtained by numerical integration, 
starting from different initial conditions. 

Analog simulation is another very convenient and known method (see [9] and Ref. therein, and [8,10,3]) 
to examine nonlinear systems. The analog device we use in this case is essentially similar to others already 
used and discussed elsewhere (see especially [9,3]). In particular, the modulation is obtained by means of 
multipliers operating in the reaction loop, in a similar way to the case of the modulated Josephson-junction 
device, discussed in detail in [3]. The values of the parameters are obtained by direct measurements on 
the experimental circuit. In particular, the damping term 5 is deduced from the resonance band width 
at the limit of small amplitudes; therefore, this measure is actually subject to some uncertainties: We 
obtain S = 0.25 with an estimated error of ±20%. The other parameters (in dimcnsionlcss units, as usual) 
are B = 1; 1/VA = 2.83; to = 1.42 ~ V2B (i.e. the frequency of the small oscillations around the 
equilibrium points of the unperturbed potential). Then, by increasing the forcing amplitude 7, we look 
for the threshold values of 7 which produce the appearance of chaotic responses, for different values of the 
amplitudes ei, £3 of the modulating perturbation. Fig. 2 shows the values we obtain for the ratio R ~ 
(where 7 are these threshold values), as a function of ei, with tz = 0. In agreement with Eq.(13), we obtain 
a completely similar behavior with ei = and varying £3. From our measures, we get the following results 



Ro = 4.3 



n = 10.1 , 



ra = 11.2 



(20) 
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to be compared with the theoretical values deduced from Eqs. (14-14') 

Ro = 2.82 , n = 11.63 , ra = 11.66 (21) 

The numerical agreement is not perfect; let us stress, however, that - apart from the uncertainties and 
unavoidable errors in the experimental determinations of the various quantities (sec [3] for some short 
comments on this point) - in any case one cannot hope that Melnikov method is able to give a precise 
determination of the threshold of chaos; rather, a common and expected result is actually that Melnikov 
theory indicates a somewhat smaller value than the threshold experimentally detected (see [6,3]). Let us 
remark on the other hand the better agreement we find for the values of ri, ra, and in particular the result 
n ~ rs (according to Eq.(14'), n = ra for w = V^), and finally the agreement shown by Fig. 2 with 
eq.(13). After these remarks, we believe that all the facts discussed up to now may be considered globally 
as a rather good test for both theory and analog experiment. 
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Figure captions 

Fig.l 

Numerical solution of Eq.(9) with different initial conditions, and A = B = u) = 1; 6 = 'y = 0.2. The 
modulation terms ei = —0.05, es = +0.04 are chosen according to the prescription (15) in order to remove 
chaos: indeed, after some short transient, the solution is periodic. 

Fig. 2 

Threshold of ehaos in the bistable potential [Eq.(9)] as a function of modulating perturbation. Here, 7 is the 
experimental value (via analog simulation) of the threshold; the ratio R = j/6 is plotted vs. the amplitude 
ei of modulation, with 63 = 0. The results with ei = and varying €3 are completely similar, see Eq.(13). 
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